Microscopic model for the magnetization plateaus in NH4CUCI3 
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A simple model consisting of three distinct dimer sublattices is proposed to describe the mag- 
netism of NH4CUCI3. It explains the occurrence of magnetization plateaus only at 1/4 and 3/4 
of the saturation magnetization. The field dependence of the excitation modes observed by ESR 
measurements is also explained by the model. The model predicts that the magnetization plateaus 
should disappear under high pressure. 

PACS numbers: 75.30.Cr, 75.10.-b, 75.40.-s 



NH 4 CuCl 3 , TICuClg, and KCuCl 3 are isostructual 
quantum spin systems consisting of two-leg ladders sep- 
arated by NHj, Tl + , and K + ions. The two-leg ladders 
are composed of Cu 2+ ions with spin 5 = 1/2 which in- 
teract antiferromagnetically (AF) through the Cl~ ions. 
These compounds can be considered as coupled two-leg 
S = 1/2 Heisenberg AF spin ladders, and show var- 
ious types of magnetization curves. Shiramura et al. 
found that NH4CUCI3 has magnetization plateaus at 1/4 
and 3/4 of the saturation moment ?i while TICUCI3 and 
KCUCI3 have no plateaus in their magnetization curves^ 
For NH4CUCI3, Kurniawan et al. performed ESR exper- 
iments which revealed that NH4CUCI3 has a finite exci- 
tation gap in the plateau regions^ By studying specific 
heat, they found AF order already at zero field and sug- 
gested a phase diagram in finite fields containing three 
magnetically ordered phases^ 

The origin of magnetization plateaus has attracted 
much interest recently. For weakly interacting dimer 
systems, the ground state is a spin singlet liquid with 
only short-range correlations between spins, and the 
triplet excitations require a finite excitation energy. An 
excited gas of triplet magnons has two characteristic 
energies^ Exchange interactions cause the transfer of 
excited triplets between neighboring dimers leading to 
a form of kinetic energy. They also lead to a short 
range repulsion between triplets in addition to the hard 
core repulsion forbidding double occupancy of a dimer. 
The triplet excitations are split in an external magnetic 
field, and the energy of the lowest component can be 
driven through zero. When the kinetic energy is dom- 
inant, which is realized in TICUCI3 and KCUCI3, we 
have no plateaus in the magnetization curved In these 
case, antiferromagnetic order perpendicular to the field 
appears simultaneously with magnetization parallel to 
the field, creating field-induced magnetic order^ which 
can be interpreted as a condensation of the lowest ly- 
ing magnon mo deA2i2ii2iiiii2ii2ii£iI£ In contrast to these 
case, when the interaction energy dominates, magne- 
tization plateaus can appear accompanied by a mag- 
netic superlattice,— Such plateaus have been observed 
at a rational fraction of the saturation momentii£ In 
SrCu2(BC>3)2, the frustrated form of the dimer lattice 
leads to a narrow bandwidth for triplet excitations and 



magnetization plateaus associated with a magnetic su- 
perlattice of localized triplets>i2iiS 

An unexpected feature of NH4CUCI3 is the absence 
of a magnetization plateau at a value 1/2 of the satu- 
ration magnetization although this would correspond to 
the simplest superlattice. In this letter, an explanation 
of the appearance of only the values 1/4 and 3/4 is pro- 
posed. Recently, Oosawa et al. performed an inelastic 
neutron scattering experiment at zero field in ND4CUCI3, 
and found two almost non-dispersive excitation branches 
at 1.8 meV and 3 meVi^ The feature of these excita- 
tion gaps is also reported in specific heat measurements 
above the transition temperature^ Since the system is 
already ordered at zero field, we can expect another low 
lying gapless branch which is hard to resolve. Very re- 
cently, Shimaoka et al. found by using NMR that 1/4 
of the dimers are in a triplet configuration below 8.5K 
in both the ordered phase and the 1/4 plateau magnetic 
field regions (H < 6T), 2J « indicating that the symmetry 
is lowered already above the magnetic transition tem- 
perature, and that three weakly interacting dimer sub- 
lattices are preformed already above the transition tem- 
perature. These results are completely different from 
the SrCu2 (1303)2 case where a phase transition breaking 
translational symmetry takes place in 1/8 plateau region. 
Motivated by these results, we propose a model composed 
of three distinct dimer sublattices to account for the mag- 
netization plateaus of NH4CUCI3. In addition, we will 
discuss the consequences of increased interladder inter- 
actions, e.g. due to application of an external pressure, 
which can lead to the suppression of the plateaus. 

The NH^~ ion is not spherical and is much larger than 
Tl + and K + ions. The NH4 molecules rotate at high tem- 
peratures, and they freeze at low temperatures S In fact, 
elastic constant shows anomaly at 70K^ and the NMR 
line splits below this temperature i2I We speculate that 
the inclusion of the NH^ ions between the ladders gives 
rise to a lattice distortion at low temperatures, modu- 
lating the exchange couplings. It is possible, because the 
angle of Cu-Cl-Cu pathway for the intradimer interaction 
is close to 90° (close to ferromagnetic exchange interac- 
tion), and a slight change of the angle may give drastic 
change in the intradimer interaction. There are various 
ways to distribute the three distinct dimer sublattices. 
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FIG. 1: Hamiltonian of our model consisting of three dis- 
tinct dimers, A, B, and C. Circles represent S = 1/2 spin. In- 
tradimer interactions are expressed by solid lines whose width 
represents the strength of the interaction. 



We show a possible simple model in Fig. 2] with three 
distinct ladders consisting of the A, B, and C dimers. 
The left and right side of each dimer is equivalent in 
this model. Two-leg ladders run along the a axis. Ji 
(i = a,b,c) is an intradimer interaction for A, B, and C 
dimers, respectively. We assume J a ~ 0.3 meV, Jb ~ 1.8 
meV, and J c ~ 3 meV to reproduce the magnon excita- 
tion modes. J[ (i — a, b, c) is an interdimer interaction 
along the ladder. Jij (i, j = a, b, c) is an interdimer inter- 
action between i and j ladders. These path ways between 
5 = 1/2 spins are extracted from the crystal structure of 
NH4CUCI3. In our model, we neglect other interactions, 
because they are expected to be small. The spin struc- 
ture in the ordered phases of NH4CUCI3 has not yet been 
determined, so we assume that its pattern is similar to 
field-induced staggered order of TICUCI3 and KCUCI3. 
Due to the modulation of the exchange couplings, the 
unit cell is larger than that of TICUCI3 and KCUCI3, 
with A and C dimers and two B dimers in the unit cell. 
In the limit A, B, and C dimers are identical, the present 
model reduces to the models for TICUCI3 and KCUCI3 in 
our previous work^ 

Taking the z axis parallel to the external magnetic 
field, we introduce the following variational wave func- 
tion at the ith dimeriSi 



i>% = c si \S) + c-fje _ix *| Tt) + ciie ixi \ jj). 



(1) 



Here, \S) is the singlet wave function. c S i, c-\i, and c^i are 
coefficients expressed as c s % — cos 0i, c-j-j = sin 0i cos fa , 
and Cw = — sin#j sin 02 . Since the left and right side of 
the dimer is equivalent in our model, the | tl) + I T4) 
triplet component does not appear in the wavefunction. 
The expectation value of the spin operator of each site of 
a dimer is given by 

(S x ) r = ~(S x )i = — -^= sin 20i (cos + sin cos 
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FIG. 2: Magnetization curves at T = OK. (a) Normalized 
staggered moment, M xy , per volume for A, B, and C dimers. 
Since there are two B dimers in the unit cell, the staggered 
magnetization for the B dimer in the phase II is about twice 
as large as for the A and C dimers in the phase I and III, 
respectively, (b) Uniform magnetization, M z , at the A, B, 
and C dimer sites, not normalized by the number of dimers. 
The total magnetization is ^4/4+_B/2+C/4, and it increases 
by 1/4, 1/2, and 1/4 upon passing through phases I, II, and 
III, respectively. Parameters were chosen to reproduce the 
experimental magnetization curves in unit of meV: J a = 0.3, 
J b = 1.75, J c = 2.95, J' a = 0.25, J' b = 0.2, J' c = 0.35, J ac = 
0.25, J bb = 0.1, Jab = 0.25, J(,c = 0.25. 



(S z ) r = (S z )i = - sin 2 



: COS 2(j)i 



(2) 



sin 20i (cos fa + sin fa ) sin x% > 



Here, (■ • • ) r m represents the expectation value on the 
right (left) side of the dimer. The perpendicular (x and 
y) component is staggered (i.e. spins are aligned oppo- 
sitely on / and r sites). The parameter x% governs the 
rotation around the z-axis, and we can determine only 
the relative phase, \a = Xc = Xb + tt, such that the 
spin configuration gains the AF interladder interaction 
energy. The angles, 9i and fa (i — a, 6, c), are variational 
parameters to be determined minimizing the expectation 
value of the Hamiltonian. This variational method is 
identical to that used in our previous paper where we 
introduced unitary transformations and minimized the c 
number term of the transformed Hamiltonian. 15 

Figure [21 shows the magnetization curves obtained by 
the variational wavefunction. There are three ordered 
phases, I, II, and III, which are driven by the A, B, and 
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FIG. 3: Schematic picture to understand the field dependence 
of excitation modes detected by ESR measurements?- wi, (J3, 
and u>6 are excitation modes from triplet to singlet, while u>2 
and u>4 are excitation modes from singlet to triplet. Wb and 
Wc represent bandwidth of triplet magnon excitations on the 
B and C dimer sublattices, respectively, which are related to 
the width of the ordered phases II and III. 
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FIG. 4: Pressure dependence of critical fields and saturation 
field. Interladder interactions (J a b, Jac, Jbb, Jbc) are increased 
by a factor of p which increases linearly with pressure: p = 
a{P — Po) + l- Here, P is pressure, Po is atmosphere pressure, 
and a is a constant. p c \ — 2.0 (p C 2 ~ 2.3) corresponds to a 
critical pressure on which the 1/4 (3/4) plateau vanishes. 



C dimers, respectively. At the A site, there is a finite 
staggered moment (M xy ) already at zero field, which in- 
duces staggered moment at B and C sites. As the field 
increases, this staggered moment first develops and then 
decreases. At H c i, which is a saturation field for the 
A dimer, the staggered moments disappear and the A 
dimer is fully polarized by the field. Since there is one 
A dimer in the unit cell as in Fig. ^ we have a magne- 
tization plateau at 1/4 for H > H c \. Above H C 2, which 
is a critical field for the B dimer, a staggered moment 
develops at the B site, inducing staggered moments at 
A and C sites. Although this field region is above the 
saturation field for the A dimer, the interaction between 
the A and B ladders induces a staggered moment at the 
A site. Accordingly, the uniform magnetization (M z ) at 
the A site decreases somewhat in the phase II as we can 
see in Fig. Hfb). Above H c s, which is a saturation field 
for the B dimer, both B and A dimers are fully polarized. 
Since there are two B dimers in the unit cell, we have a 
magnetization plateau at 3/4 for H > H C 3. Above H C 4 
(critical field for the C dimer), a staggered moment de- 
velops at the C site, inducing the staggered moment at A 
and B sites. Above H s (saturation field for the C dimer), 
the all dimers are fully polarized. Thus, the magnetiza- 
tion plateaus at 1/4 and 3/4 of the saturation moment 
can be understood as successive quantum phase transi- 
tions driven by magnetic field, which occur in interacting 
three distinct dimer sublattices. 

Specific heat and high-field magnetization measure- 
ments found three magnetically ordered phases on the 
H — T phase diagram^ It is interesting to compare it 
with Fig. Ufa). Since we expect a larger staggered mag- 
netic moment leads to a higher transition temperature, 
we can account for the resemblance between the H — T 
phase diagram and M xy (H) curves. 



In our model, the excitation branches at 1.8 meV and 
3 meV found by neutron scattering measurements are 
triplet excitations. Therefore, each of them should split 
into three branches when we apply external magnetic 
field. It can explain the origins of the four low lying 
excitation modes, (i = 1 ~ 4) (see Fig. |3J|, detected 
by ESR measurements^ In the 1/4-plateau region, a w\- 
mode increases linearly with the field, while the a>2-mode 
decreases and goes soft at H c2 . In this field region, each 
A dimer is fully polarized by the field, and is in the lowest 
lying triplet | ||), while the configuration of B dimer is 
dominated by the singlet component. Therefore, the uj\- 
mode can be identified with an excitation from a triplet 
to singlet on the A dimer, and the W2-mode with an ex- 
citation from the singlet to the lowest lying triplet on 
the B dimer. Thus, an excitation energy gap opens up 
in the 1/4-plateau region. In the 3/4-plateau region, the 
energy of the W3-mode increases with the field, while U4- 
mode decreases and goes soft at H C 4. Similarly to the 
1/4-plateau region, the W3-mode can be identified with 
an excitation from the lowest triplet to the singlet on the 
B dimer, and L04 with an excitation from singlet to the 
lowest triplet on the C dimer. luq in Fig. [21 is an excita- 
tion from the lowest triplet to the singlet on the C dimer, 
which is not detected by ESR measurements. There is 
no singlet-triplet excitation gap on the A dimer, since an 
AF order appears on the A dimer already at zero field. 

Consider now what we can predict from our model. 
The critical fields (H c i, H C 2, H C 3, H c a) are functions 
of interdimer interactions which we can control by pres- 
sure. If we assume that interladder interactions increase 
linearly with pressure, while the intraladder interactions 
are constant, we obtain the pressure dependence of the 
critical fields shown in Fig. 01 As pressure increases, 
the three distinct dimer sublattices couple more strongly, 
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FIG. 5: Magnetization curves under pressure p — p c \. (a) 
Staggered moment, (b) uniform magnetization. Interladder 
interactions (J ac , Jbb, Jab, Jbc) are simply increased by a fac- 
tor of pd. The other couplings are the same as in Fig. [5] 

and a stronger field is required to saturate spins. Thus 
the plateau onset fields, H c t, H c ^, and H s , increase with 
pressure, since they correspond to saturation fields for 
A, B, and C dimers, respectively. The plateau end fields, 
H C 2 and H C 4, are critical fields where the lowest triplet 
component on the B and C dimers are driven to zero 
energy. Consider first the case of the B-sublattice. The 
increase with pressure of the intrasublattice, Jbb, broad- 
ens the triplet magnon bandwidth for the B dimer and 
lowers H C 2- The coupling through the C-sublattice Jbc 
acts similarly but the coupling to the already polarized 
A-sublattice acts oppositely to increase H c2 , since the 
spins on the A sites are oriented parallel to the field. The 
net result however is a decrease in H C 2(p) as in Fig. 0] 
This decrease leads to the disappearance of the plateau 
between H c \ and H C 2 at a critical pressure. In Fig. 
we show the magnetizations M xy (H) and M Z (H) at this 
quantum critical point. The case of the 3/4-plateau is 
slightly different, since H c i(p) increases with the pres- 
sure. This occurs because there is negligible direct in- 
trasublattice coupling for C dimers and both the A and 
B sublattices are now polarized parallel to H. None the 
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less the 3/4-plateau width also narrows with increasing 
p and eventually disappears at a higher critical pressure 
above which all the ordered phases I, II, and III merge 
into a single phase. 

In summary, we have proposed a model consisting of 
three distinct dimer sublattices in order to account for the 
magnetization plateaus that appear in NH4CUCI3. There 
are various way to distribute the three distinct dimer sub- 
lattices. In this paper, we proposed a simple model in 
which the dimer sublattices distribute on the b — c plane 
as illustrated in Fig. ^ In this case, triplet magnon ex- 
citations propagate predominantly along the ladder di- 
rection. Below H C 2, the spin structure of our model is 
characterized by a wave vector (0,7r,7r) or (0,0,27r). An 
alternate possibility is to distribute the three distinct 
dimer sublattices along the ladders (a axis), which is pro- 
posed by studying NMRiSi Our theory is applicable also 
to this case. The distribution pattern of the three dis- 
tinct dimer sublattices is not essential to our main results. 
The important assumption of our model is that there are 
three weakly interacting distinct dimer sublattices whose 
intradimer interactions are characterized by the values 
0.3 meV, 1.8 meV, and 3 meV observed by the neutron 
experiments^, The volume fraction of the dimers should 
be 1/4, 1/2, and 1/4, respectively. Note this is consistent 
with the fact that the intensity of the magnon branch 
at 1.8 meV is about twice as large as the branch at 3 
meV)«i and it explains why the magnetization plateaus 
take place only at 1/4 and 3/4. We made a prediction 
that the magnetization plateaus disappear by applying 
pressure, which is supported by a recent magnetization 
measurement under high pressured 
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